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>- ■ 1 Introduction 

\o 

■ Magnetic fields permeate the Universe. They are found in planets, stars, accre- 
I tion discs, galaxies, clusters of galaxies, and the intergalactic medium. While 

■ there is often a component of the field that is spatially coherent at the scale 
I of the astrophysical object, the field lines are tangled chaotically and there 

■ are magnetic fluctuations at scales that range over orders of magnitude. The 
fH I cause of this disorder is the turbulent state of the plasma in these systems. 
pi i' This plasma is, as a rule, highly conducting, so the magnetic field lines are en- 
^ I trained by (frozen into) the fluid motion. As the fields are stretched and bent 
^ • by the turbulence, they can resist deformation by exerting the Lorentz force 

I on the plasma. The turbulent advection of the magnetic field and the field's 

■ back reaction together give rise to the statistically steady state of fully devel- 
oped MHD turbulence. In this state, energy and momentum injected at large 
(object-size) scales are transfered to smaller scales and eventually dissipated. 

. , Despite over fifty years of research and many major advances, a satisfactory 

■ theory of MHD turbulence remains elusive. Indeed, even the simplest (most 
idealised) cases are still not fully understood. One would hope that there are 
universal properties of MHD turbulence that hold in all applications — or 
at least in a class of applications. Among the most important questions for 
astrophysics that a successful theory of turbulence must answer are: 

• How does the turbulence amplify, sustain and shape magnetic fields? What 
is the structure and spectrum of this field at large and small scales? The 
problem of turbulence in astrophysics is thus directly related to the fun- 
damental problem of magnetogenesis. 
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• How is energy cascaded and dissipated in plasma turbulence? In accretion 
discs and the solar corona, for example, one would like to know if the 
turbulence heats ions or electrons predominantly |52j . 

• How does the turbulent flow and magnetic field enhance or inhibit the 
transport of heat, (angular) momentum, and cosmic rays? Again in accre- 
tion discs, a key parameter is the effective turbulent viscosity that causes 
the transport of angular momentum and enables accretion [62] ■ In cluster 
physics, an understanding of how viscous heating and thermal conduction 
in a turbulent magnetised plasma balance the radiative cooling is necessary 
to explain the observed global temperature profiles [12] . 

In this chapter, we discuss the current understanding of the most basic 
properties of astrophysical MHD turbulence. We emphasise possible universal 
aspects of the theory. We shall touch primarily on two applications: turbu- 
lence in the solar wind and in clusters of galaxies. These are, in a certain (very 
approximate) sense, two "pure" cases of small-scale turbulence, where theo- 
retical models of the two main regimes of MHD turbulence (discussed in §[2] 
and §[3]) can be put to the test. They are also good examples of a complication 
that is more or less generic in astrophysical plasmas: the MHD description is, 
in fact, insufficient for astrophysical turbulence and plasma physics must make 
an entrance. Why this is so will be explained in §|4j 

The astrophysical plasma turbulence is even more of a terra incognita than 
the MHD turbulence, so we shall start with the equations of incompressible 
MHD — the simplest equations that describe (subsonic) turbulent dynamics 
in a conducting medium: 

-Vp + uAu + B-VB + i, V-u = 0, (1) 

B-Vu-|-r;AB, (2) 

where u is the velocity field, d/dt ^ dt + u-'V the convective derivative, p the 
pressure (scaled by the constant density p and determined by the incompress- 
ibility constraint), B the magnetic field scaled by {Airpy^^, v the kinematic 
viscosity, r] the magnetic diffusivity, and f the body force that models large- 
scale energy input. The specific energy injection mechanisms vary: typically, 
in astrophysics, these are either background gradients (e.g., the temperature 
gradient in stellar convective zones, the Keplerian velocity shear in accretion 
discs), which mediate the conversion of gravitational energy into kinetic en- 
ergy of fluid motion, or direct sources of energy such as the supernovae in the 
interstellar medium or active galactic nuclei in galaxy clusters. What all these 
injection mechanisms have in common is that the scale at which they oper- 
ate, hereafter denoted by L, is large, comparable with the size of the system. 
While the large-scale dynamics depend on the specific astrophysical situation, 
it is common to assume that, once the energy has cascaded down to scales 
substantially smaller than L, the nonlinear dynamics are universal. 
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The universality of small scales is a cornerstone of all theories of turbu- 
lence. It goes back to Kolmogorov's 1941 dimensional theory (or K41 [35]; see 
[37], §33 for a lucid and concise exposition). Here is an outline Kolmogorov's 
reasoning. Consider Eq. ([1]) without the magnetic term. Denote the typical 
fluctuating velocity difference across scale L by 5ul. The energy associated 
with these fluctuations is 5u\ and the characteristic time for this energy to cas- 
cade to smaller scales by nonlinear coupling is L/5ul- The total specific power 
(energy flux) going into the turbulent cascade is then e = (u • f) 5u\/L. 
In a statistically stationary situation, all this power must be dissipated, so 
e = z^(|Vup). Since e is a finite quantity completely defined by the large-scale 
energy-injection process, it cannot depend on v. For very small z/, this im- 
plies that the velocity must develop very small scales so that z/(|Vup) has 
a constant limit as v ^ +0. The only quantity with dimensions of length 
that one can construct out of e and v is ^ (y^ jt)^!'^ ^ Re~'^^^L, where 
Re ~ 6ulL/i> is the Reynolds number. In astrophysical applications. Re is 
usually large, so the viscous dissipation occurs at scales li, ^ L. The energy 
injected at the large scale L must be transfered to the small scale li, across 
a range of scales (the inertial range). The hydrodynamic turbulence theory 
assumes that the physics in this range is universal, i.e., it depends neither 
on the energy-injection mechanism nor on the dissipation mechanism. Four 
further assumptions are made about the inertial range: homogeneity (no 
special points), scale invariance (no special scales), isotropy (no special 
directions), and locality of interactions (interactions between comparable 
scales dominate). Then, at each scale / such that L ^ I ^ l,^, the total power 
e must arrive from larger scales and be passed on to smaller scales: 

er^Suf/n, (3) 

where 6ui is the velocity difference across scale I and t; the cascade time. 
Dimensionally, only one time scale can be constructed out of the local quan- 
tities Sui and I: ti ^ l/Sui. Substituting this into Eq. ([3]) and solving for Sui, 
we arrive at Kolmogorov's scaling: 6ui ~ {el^^^, or, for the energy spectrum 
Eik), 

Suf ^ / dk'Eik') ^ £2/3^-2/3 ^ ^(fc) _ e2/3fc-5/3_ (4) 

Jk=l/l 

The history of the theory of MHD turbulence over the last half century has 
been that of a succession of attempts to adapt the K41-style thinking to fluids 
carrying magnetic fields. In the next two sections, we give an overview of these 
efforts and of the resulting gradual realisation that the key assumptions of 
the small-scale universality, isotropy and locality of interactions fail in various 
MHD contexts. 
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2 Alfvenic turbulence 

Let us consider the case of a plasma threaded by a straight uniform magnetic 
field Bq of some external (i.e., large-scale) origin. Let us also consider weak 
forcing so that the fundamental turbulent excitations are small-amplitude 
wave-like disturbances propagating along the mean field. We will refer to such 
a limit as Alfvenic turbulence — it is manifestly anisotropic. 

2.1 Iroshnikov Kraichnan turbulence 

If we split the magnetic field into the mean and fluctuating parts, B = Bo-l-(5B, 
and introduce Elsasser [Mj variables = u ± (5B, Eqs. ([T]) and ^ take a 
symmetric form: 

atz± T waV||Z± + z^ • Vz± = -Vp + ^^-i-^Az± + ^^y^Az^ + f , (5) 

where wa = |Bo| is the Alfven speed and V|| is the gradient in the direction of 
the mean field Bq. The Elsasser equations have a simple exact solution: if z+ = 
or z^ = 0, the nonlinear term vanishes and the other, non-zero, Elsasser field 
is simply a fluctuation of arbitrary shape and magnitude propagating along 
the mean field at the Alfven speed va- Kraichnan [Mj realised in 1965 that 
the Elsasser form of the MHD equations only allows nonlinear interactions 
between countcrpropagating such fluctuations. The phcnomcnological theory 
that he and, independently, Iroshnikov [30| . developed on the basis of this 
idea (the IK theory) can be summarised as follows. 

Following the general philosophy of K41, assume that only fluctuations 
of comparable scales interact (locality of interactions) and consider these 
interactions in the inertial range, comprising scales I smaller than the forcing 
scale L and larger than the (still to be determined) dissipation scale. Let us 
think of the fluctuations propagating in either direction as trains of spatially 
localised Alfven- wav^ packets of parallel (to the mean field) extent ^|| and 
perpendicular extent I (we shall not, for the time being, specify how /y relates 
to I). Assume further that Jz;"^ ^ Szf ~ 6ui ^ SBi. We can again use Eq. 
([31) for the energy flux through scale /, but there is, unlike in the case of 
purely hydrodynamic turbulence, no longer a dimensional inevitability about 
the determination of the cascade time r/ because two physical time scales 
are associated with each wave packet: the Alfven time ta{1) ^ Iw/va and the 
strain (or "eddy") timeTs(^) l/Sui. To state this complication in a somewhat 
more formal way, there are three dimcnsionlcss combinations in the problem 
of MHD turbulence: el/Suf, Sui/va and so the dimensional analysis does 
not uniquely determine scalings and further physics input is needed. 

^ Waves in incompressible MHD can have either the Alfven- or the slow-wave polar- 
isation. Since both propagate at the Alfven speed, we shall, for simplicity, refer 
to them as Alfven waves. The differences between the Alfven- and slow-wave 
cascades are explained in detail at the end of §[2131 



Turbulence and Magnetic Fields in Astrophysical Plasmas 5 



Two countcrpropagating wave packets take an Alfvcn time to pass through 
each other. During this time, the amphtude of either packet is changed by 

A5ui ^ ^^TA ^ 6ui — . (6) 

The IK theory now assumes weak interactions, ASui <C 6ui ta ^ Tg. 
The cascade time ti is estimated as the time it takes (after many interactions) 
to change Sui by an amount comparable to itself. If the changes in amplitude 
accumulate like a random walk, we have 

t I — ; — o 72 

E.J. - Ta / I (. r , T I, VA 

Adul ~ dui — \ — dui tor t n ti ^ — , ^ 9 ■ (' ) 

Ts V TA TA /jjO-uf 

Substituting the latter formula into Eq. ([3]), we get 

5ui ^ {evA)"\"'l"^. (8) 

The final IK assumption, which at the time seemed reasonable in light of the 
success of the K41 theory, was that of isotropy, fixing the dimensionlcss ratio 
~ 1, and, therefore, the scaling: 

5ui ^ {evAf'^l^'^ E{k) ^ (ewA)'/'fc-'/'. (9) 



2.2 Turbulence in the solar wind 

The solar wind, famously predicted by Parker [SD], was the first astrophysical 
plasma in which direct measurements of turbulence became possible [11) . A 
host of subsequent observations (for a concise review, see [24] ) revealed pow- 
erlike spectra of velocity and magnetic fluctuations in what is believed to be 
the inertial range of scales extending roughly from 10^ to 10'^ km. The mean 
magnetic field is Bq ^ 10. ..10^ /iG, while the fluctuating part SB is a factor 
of a few smaller. The velocity dispersion is 6u ~ 10^ km/s, approximately in 
energy equipartition with SB. The u and B fluctuations are highly correlated 
at all scales and almost undoubtedly Alfvenic. It is, therefore, natural to think 
of the solar wind as a space laboratory conveniently at our disposal to test 
theories of Alfvenic turbulence in astrophysical conditions. 

For nearly 30 years following Kraichnan's paper [53], the IK theory was 
accepted as the correct extension of K41 to MHD turbulence and, therefore, 
with minor modifications allowing for the observed imbalance between the 
energies of the z+ and z fluctuations [13], also to the turbulence in the 
solar wind. However, alarm bells were sounding already in 1970s and 1980s 
when measurements of the solar-wind turbulence revealed that it was strongly 
anisotropic with Z|| > l± [5] (see also [33]) and that its spectral index was 
closer to -5/3 than to -3/2 [t^FI Numerical simulations have confirmed the 
anisotropy of MHD turbulence in the presence of a strong mean field [40l[9l[46]. 

Another classic example of a —5/3 scaling in astrophysical turbulence is the 
spectrum of electron density fluctuations (thought to trace the velocity spectrum) 
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2.3 Weak turbulence 



The realisation that the isotropy assumption must be abandoned led to a 
reexamination of the Alfven-wave interactions in MHD turbulence. If the as- 
sumption of weak interactions is kept, MHD turbulence can be regarded as 
an ensemble of waves, whose wavevectors k and frequencies uj^{\i) = ±A:||Ua 
have to satisfy resonance conditions in order for an interaction to occur. For 
three-wave interactions (1 and 2 counterpropagating, giving rise to 3), 

ki+k2=k3 => +fc||2 = fc||3, (10) 

tJ±(ki) +CjT(k2) =C^±(k3) fc||i-fc||2 = fc||3, (11) 

whence fc||2 — and fc||3 — fc||i. Thus (i) interactions do not change fc||; (ii) 
interactions arc mediated by modes with fc|| = 0, which arc quasi-2D fluctua- 
tions rather than waves [45l [63l H 

The first of these conclusions suggests a quick fix of the IK theory: take 
/|| ^ fc||p^ — const (the wavenumber at which the waves are launched) and 
/ ^ lj_ in Eq. ^ (no parallel cascade). Then the spectrum is [55] 

E{k^) ^ (efc||oi'A)'/'fcl'. (12) 

The same result can be obtained via a formal calculation based on the stan- 
dard weak-turbulence theory [20l[39j. However, it is not uniformly valid at all 
k±. Indeed, let us check if the assumption of weak interactions, ta Ts, is 
actually satisfied by the scaling relation ([5]) with /|| ^ k^^: 



— ~ ^TiTT^ « 1 ^ » ^* = — ~ -?rijrT^ (13) 



2 r ' 



where 6ul is the velocity at the outer scale (the rms velocity). Thus, if Re = 
SulL/v and Rm = SulL/tj are large enough, the inertial range will always 
contain a scale h below which the interactions arc no longer weak|f| 



in the interstellar medium — the famous "power law in the sky," which appears 
to hold across 12 decades of scales [T]. 

^ Goldreich & Sridhar [23) argued that the time it takes three waves to realise that 
one of them has zero frequency is infinite and, therefore, the weak-interaction 
approximation cannot be used. This difficulty can, in fact, be removed by noticing 
that the fcy = modes have a finite correlation time, but we do not have space 
to discuss this rather subtle issue here (two relevant references are [20ll39j ). 

^ There is also an upper limit to the scales at which Eq. (|12p is applicable. The 
boundary conditions at the ends of the "box" are unimportant only if the cascade 
time (0 is shorter than the time it takes an Alfvenic fluctuation to cross the box: 
Ti <c; L\\/va ^ L* = (£/fe||o«A)^''^-^|| , where is the length of the box 

along the mean field. Demanding that L» > L, the perpendicular size of the box, 
we get a lower limit on the aspect ratio of the box: > k\\(^L[vA/ 5ulY . If this 
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R. S. Iroshnikov (1937-1991) R. H. Kraichnan 




P. Goldreich S. Sridhar 



Fig. 1. IK and GS (photo of R. S. Iroshnikov courtesy of Sternberg Astronomical 
Institute; photo of R. H. Kraichnan courtesy of the Johns Hopkins University). 

2.4 Goldreich Sridhar turbulence 

In 1995, Goldreich & Sridhar [22] conjectured that the strong turbulence below 
the scale U should satisfy 

TA ~ T"s <^ va/5ui, (14) 

a property that has come to be known as the critical balance. Goldreich & 
Sridhar argued that when ta ^ Ts, the weak turbulence theory "pushes" the 
spectrum towards the approximate equality (jl4p . They also argued that when 
Ta S> Tg, motions along the field lines are decorrelated and naturally develop 
the critical balance. 

The critical balance fixes the relation between two of the three dimension- 
less combinations in MHD turbulence. Since now there is only one natural time 
scale associated with fluctuations at scale I, this time scale is now assumed 
to be the cascade time, r; ^ Tg. This brings back Kolmogorov's spectrum ^ 
for the perpendicular cascade. The parallel cascade is now also present but is 
weaker: from Eq. (fT4|) . 

^ll~^-A6-l/3/f (15) 

The scalings (|4|) and should hold at all scales l± <C h and above 
the dissipation scale: either viscous l,^ ^ (z^^/e)^/"* or resistive Iri ~ (ry'^/e)^/^, 

is not satisfied, the physical (non-periodic) boundary conditions may impose a 
limit on the perpendicular field-line wander and thus effectively forbid the fcy — 
modes. It has been suggested [23] that a weak-turbulence theory based on 4-wave 
interactions [B^ should then be used at l± > L,. 
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whichever is larger. Comparing these scales with [Eq. (fT5|) ]. we note that the 

strong-turbulence range is non-empty only if Re, Rm ^ (^k^^oL)^ (va/Sul)^ , a 
condition that is effortlessly satisfied in most astrophysical cases but should 
be kept in mind when numerical simulations are undertaken. 

The Goldreich-Sridhar (GS) theory has now replaced the IK theory as the 
standard accepted description of MHD turbulence. The feeling that the GS 
theory is the right one. created by the solar wind [42^ and ISM [1] observations 
that show a k~^^^ spectrum, is, however, somewhat spoiled by the consistent 
failure of the numerical simulations to produce such a spectrum [401 146] . In- 
stead, a spectral index closer to IK's —3/2 is obtained (this seems to be the 
more pronounced the stronger the mean field), although the turbulence is 
definitely anisotropic and the GS relation appears to be satisfied [101 [3] ! 
This trouble has been blamed on intermittency [40j , a perennial scapegoat of 
turbulence theory, but a non-speculative solution remains to be found. 

The puzzling refusal of the numerical MHD turbulence to agree with ei- 
ther the GS theory or, indeed, with the solar-wind observations highlights 
the rather shaky quality of the existing physical understanding of what really 
happens in a turbulent magnetic fluid on the dynamical level. One conceptual 
difference between MHD and hydrodynamic turbulence is the possibility of 
long-time correlations. In the large- i?™ limit, the magnetic field is determined 
by the displacement of the plasma, i.e., the time integral of the (Lagrangian) 
velocity. In a stable plasma, the field-line tension tries to return the field line 
to the unperturbed equilibrium position. Only "interchange" (A:|| = 0) mo- 
tions of the entire field lines are not subject to this "spring-back" effect. Such 
motions are often ruled out by geometry or boundary conditions (cf. footnote 
[6]). Thus, fluid elements in MHD cannot simply random walk as this would 
increase (without bound) the field-line tension. However, they may random 
walk for a substantial period before the tension returns them back to the 
equilibrium state. The role of such long-time correlations in MHD turbulence 
is unknown. 



Reduced MHD, the decoupling of the Alfven-wave cascade, and turbu- 
lence in the interstellar medium. We now give a rigourous demonstration of 
how the turbulent cascade associated with the Alfven waves (or, more precisely, 
Alfven-wave-polarised fluctuations) decouples from the cascades of the slow waves 
and entropy fluctuations. Let us start with the equations of compressible MHD: 

f^-pV.u, (16) 

p5^-vfp+f)+^, (17) 
dt \ Stt y 47r 

ds p 5 , , 

— = 0, S=— , 7=TT, (18 

^=B-Vu-BV-u. (19) 
dt ^ ' 
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Consider a uniform static equilibrium with a straight magnetic field, so p = po + 5p, 
p = po + 5p, B = Bo + 5B. Based on observational and numerical evidence, it 
is safe to assume that the turbulence in such a system will be anisotropic with 
fcy <Si k±. Let us, therefore, introduce a small parameter e ~ k^\/k± and carry out 
a systematic expansion of Eqs. (|16|l - (|19| l in e. In this expansion, the fluctuations 
are treated as small, but not arbitrarily so: in order to estimate their size, we shall 
adopt the critical-balance conjecture HJ), which is now treated not as a detailed 
scaling prescription but as an ordering assumption. This allows us to introduce the 
following ordering: 

^^ILL^!iL^^^ ^ ^ ^ e (20) 

Po VA VA po Bo Bo k± 

where we have also assumed that the velocity and magnetic-fleld fluctuations have 
the character of Alfven and slow waves (JB ~ u) and that the relative amplitudes of 
the Alfven- wave-polarised fluctuations {u±/va, SB±/Bo), slow- wave-polarised fluc- 
tuations 5B\\/Bo) and density fluctuations (Sp/po) are all the same orderQ 
We further assume that the characteristic frequency of the fluctuations is ~ fcy^A, 
which means that the fast waves, for which cj ~ k±^v\ -f ci , where Cs — 7Po/po is 
the sound speed, are ordered out. 

We start by observing that the Alfven-wave-polarised fluctuations are two- 
dimensionally solenoidal: since V ■ u = 0(6^) [from Eq. ^] and V • B = 0, 
separating the 0(e) part of these divergences gives Vx ■ ux = Vx ■ S'B± = 0. We 
may, therefore, express ux and 5Bx in terms of scalar stream (flux) functions: 

ux = bo X Vx0, -7=^ = bo X Vx^, (21) 
v47rpo 

where bo = Bo /-Bo- Evolution equations for and tp are obtained by substituting 
the expressions (|2ip into the perpendicular parts of the induction equation (|19[) and 
the momentum equation (|17p — of the latter the curl is taken to annihilate the 
pressure term. Keeping only the terms of the lowest order, O(e^), we get 

+ {<l,,^} = va\7 ii<P, (22) 

^Vi0 + {0, Vi0} = i;aV|| Vi^/> + {i>, Vli^} , (23) 

where {(j},tp} — bo ■ (V ±(t> x Vx'/') and, to lowest order, 

^ = |+ux.Vx = |+{0,..-}, (24) 

g.V^V„ + ^.Vx^V„ + ^W,...}. (25) 

Eqs. (|22|) and (|23|l are known as the Reduced Magnetohydrodynamics (RMHD). 
They were flrst derived by Strauss [66] in the context of fusion plasmas. They form 



^ Strictly speaking, whether this is the case depends on the energy sources that 
drive the turbulence: as we are about to see, if no slow waves are launched, none 
will be present. However, it is safe to assume in astrophysical contexts that the 
large-scale energy input is random and, therefore, comparable power is injected 
in all types of fluctuations. 
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a closed set, meaning that the Allven-wave cascade decouples from the slow waves 
and density fluctuations. 

In order to derive evolution equations for the latter, let us revisit the perpendic- 
ular part of the momentum equation and use Eq. (|20p to order terms in it. In the 
lowest order, 0(e), we get the pressure balance 

4tt } po ci Bo 



Using Eq. p6|l and the entropy equation (|18|) . we get 

^^^0, ^ = ^-^^ = -yf^ + 4^y (27) 
at So So Po po V po c| Bo J 

where so ~ Po/Po- On the other hand, from the continuity equation (|16p and the 
parallel component of the induction equation (|19p , 

dt \po Bo J Bo 



Combining Eqs. p7|) and pSf) . we obtain 
d Sp 1 



dtpo l + ci/«iBo ^ ^ 

dSB« IB 

Finally, we take the parallel component of the momentum equation (|17|l and notice 
that, due to Eq. (|26|l and to the smallness of the parallel gradients, the pressure 
term is 0(e'^), while the inertial and tension terms are O(e^). Therefore, 

Eqs. pO|l and (|3H) describe the slow-wave-polarised fluctuations, while Eq. (|27|) de- 
scribes the zero-frequency entropy mode. The nonlinearity in these equations enters 
via the derivatives defined in Eqs. (|24|l - (|25|l and is due solely to interactions with 
Alfven waves. 

Naturally, the reduced equations derived above can be cast in the Elsasser form. 
If we introduce Elsasser potentials (^^ = (f) zL ip, Eqs. (|22|l and (|23|l become 

|vic^ T^^V||Vic± = -i [{c+, vir} + {r, vic+} T vi {c+,r}] • (32) 

This is the same as the perpendicular part of Eq. |[5]) with = bo x Vx(^*. The 
key property that only counterpropagating Alfven waves interact is manifest here. 
For the slow-wave variables, we may introduce generalised Elsasser fields; 



SB« 



V47rpo V 



2 \ 1/2 

1 + ^ ■ (33) 



Straightforwardly, the evolution equation for these fields is 
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This equation reduces to the parallel part ol Eq. © in the limit va ^ Cs- This 
is known as the high-/? limit, with the plasma beta defined hy P — Stvpo/Bq — 
(2/7)Cs/tj^. We see that only in this limit do the slow waves interact exclusively 
with the counterpropagating Alfven waves. For general /3, the phase speed of the 
slow waves is smaller than that of the Alfven waves and, therefore, Alfven waves 
can "catch up" and interact with the slow waves that travel in the same direction. 

In astrophysical turbulence, f3 tends to be moderately high[f| for example, in the 
interstellar medium, /3 ~ 10 by order of magnitude. In the high-/3 limit, which is 
equivalent to the incompressible approximation for the slow waves, density fluctua- 
tions are due solely to the entropy mode. They decouple from the slow-wave cascade 
and are passively mixed by the Alfven- wave turbulence: dSp/dt = [Eq. (|27|l or 
Eq. (|29|) . Cs ^ Va]- By a dimensional argument similar to K41, the spectrum of 
such a field is expected to follow the spectrum of the underlying turbulence [IS]: in 
the GS theory, k^^^^. It is precisely the electron-density spectrum (deduced from 
observations of the scintillation of radio sources due to the scattering of radio waves 
by the interstellar medium) that provides the evidence of the k'^''^^ scaling in the 
interstellar turbulence [T] . The explanation of this density spectrum in terms of pas- 
sive mixing of the entropy mode, originally conjectured by Higdon [23, is developed 
on the basis of the GS theory in Ref. [38] . 

Thus, the anisotropy and critical balance taken as ordering assumptions lead to a 
neat decomposition of the MHD turbulent cascade into a decoupled Alfven- wave cas- 
cade and cascades of slow waves and entropy fluctuations passively scattered/mixed 
by the Alfven waveslf] The validity of this decomposition and, especially, of the 
RMHD equations (|22|l and (|23|) turns out to extend to collisionless scales, where the 
MHD equations (|16|) - (|19p cannot be used: this will be briefly discussed in ij |4.3l 



3 Isotropic MHD turbulence 

Let us now consider the case of isotropic MHD turbulence, i.e., a turbulent 
plasma where no mean field is imposed externally. Can the scaling theories 
reviewed above be adapted to this case? A popular view, due originally to 
Kraichnan [34] , is that everything remains the same with the role of the mean 
field Bq now played by the magnetic fluctuations at the outer scale, 5i?L, 



The solar corona, where j3 ~ 10^ , is one prominent exception. 
Eqs. p7|) . (|32|) and (|34p imply that, at arbitrary /3, there are five conserved quan- 
tities: la = {|5s|^) (entropy fluctuations), I"^ = {|V("*p) (right /left-propagating 
Alfven waves), — (right/left-propagating slow waves). and are 

always cascaded by interaction with each other. Is is passively mixed by and 
/J, 7|i'^ are passively scattered by and, unless /3 3> 1, also by ij^. 
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Fig. 2. Cross sections ol the absolute values of u (left panel) and B (right panel) in 
the saturated state of a simulation with Re ~ 100, Prm = 10 (run B of Ref. [57] )■ 



while at smaller seales, the turbulence is again a local (in scale) cascade of 
Alfvenic {5ui ^ 5Bi) fluctuations. This picture is only plausible if the mag- 
netic energy is dominated by the outer-scale fluctuations, an assumption that 
does not appear to hold in the numerical simulations of forced isotropic MHD 
turbulence [S7]. Instead, the magnetic energy is concentrated at small scales, 
where the magnetic fluctuations signiflcantly exceed the velocity fluctuations, 
with no sign of the scale-by-scale equipartition implied for an Alfvenic cas- 
cade0 These features are especially pronounced when the magnetic Prandtl 
number Pr-^ = u/r] = Rm/Re ^ 1, i.e., when the magnetic cutoff scale lies 
below the viscous cutoff of the velocity fluctuations (Fig. [2]). The numerically 
more accessible case of Pr,„ > 1, while non-asymptotic and, therefore, harder 
to interpret, retains most of the features of the large- Pr^ regime. A handy 
formula for Pr^ based on the Spitzer [64j values of v and 77 for fully ionised 
plasmas is 

Prra - lO-^rV", (35) 

where T is the temperature in Kelvin and n the particle density in cm^'^. Eq. 
([35]) tends to give very large values for hot diffuse astrophysical plasmas: e.g., 
10^^ for the warm interstellar medium, 10^^ for galaxy clusters. 

Let us examine the situation in more detail. In the absence of a mean 
field, all magnetic fields are generated and maintained by the turbulence itself, 

^° This is true for the case of forced turbulence. Simulations of the decaying case [4] 
present a rather different picture: there is still no scale-by-scale equipartition but 
the magnetic energy heavily dominates at the large scales — most likely due to a 
large-scale force-free component controlling the decay. The difference between the 
numerical results on the decaying and forced MHD turbulence points to another 
break down in universality in stark contrast with the basic similarity of the two 
regimes in the hydrodynamic case. 
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Fig. 3. Stretching/shearing a magnetic-field line. 



i.e., isotropic MHD turbulence is the saturated state of the turbulent (small- 
scale) dynamo. Therefore, we start by considering how a weak (dynamically 
unimportant) magnetic field is amplified by turbulence in a large- Pr,„ MHD 
fluid and what kind of field can be produced this way. 

3.1 Small-scale dynamo 

Many specific deterministic flows have been studied numerically and analy- 
tically and shown to be dynamos [5] . While rigourously determining whether 
any given flow is a dynamo is virtually always a formidable mathematical 
challenge, the combination of numerical and analytical experience of the last 
fifty years suggests that smooth three-dimensional flows with chaotic trajec- 
tories tend to have the dynamo property provided the magnetic Reynolds 
number exceeds a certain threshold, Rm > Rm.c ~ 10^. ..10^. In particular, 
the ability of Kolmogorov turbulence to amplify magnetic fields is a solid nu- 
merical fact first established by Meneguzzi et al. in 1981 |44] and since then 
confirmed in many numerical studies with ever-increasing resolutions (most 
recently [57l[27]). It was, in fact, Batchelor who realised already in 1950 [2] 
that the growth of magnetic fluctuations in a random fiow should occur sim- 
ply as a consequence of the random stretching of the field lines and that it 
should proceed at the rate of strain associated with the fiow. In Kolmogorov 
turbulence, the largest rate of strain ~ Sui/l is associated with the smallest 
scale I ^ Iv — the viscous scale, so it is the viscous-scale motions that dom- 
inantly amplify the field (at large Pr^). Note that the velocity field at the 
viscous scale is random but smooth, so the small-scale dynamo in Kolmogorov 
turbulence belongs to the same class as fast dynamos in smooth single-scale 



A repeated application of random stretching/shearing to a tangled mag- 
netic field produces direction reversals at arbitrarily small scales, giving rise 
to a folded field structure (Fig. [3]). It is an essential property of this structure 
that the strength of the field and the curvature of the field lines are anticor- 
related: wherever the field is growing it is relatively straight (i.e., curved only 
on the scale of the fiow), whereas in the bending regions, where the curvature 
is large, the field is weak. A quantitative theory of the folded structure can 



fiows [i[57]. 
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Fig. 4. Ya. B. Zeldovich (1914-1987) (photo courtesy of M. Ya. Ovchinnikova) . 



be constructed based on the joint statistics of the field strength B = |B| and 
curvature K = b • Vb, where b = B/i? [SS]0 The curvature is a quantity 
easily measured in numerical simulations, which confirm the overall straight- 
ness of the field and the curvature-field-strength anticorrelation [ST]- At the 
end of this section, we shall give a simple demonstration of the validity of the 
folded structure. 

The scale of the direction reversals is limited from below only by Ohmic 
diffusion: for Kolmogorov turbulence, balancing the rate of strain at the vis- 
cous scale with diffusion and taking Pr^ ^ 1 gives the resistive cutoff 

5uiJl,^iT/ll ^ l,,^Pr;^^'X. (36) 

If random stretching gives rise to magnetic fields with reversals at the re- 
sistive scale, why are these fields not eliminated by diffusion? In other words, 
how is the small-scale dynamo possible? It turns out that, in 3D, there are 
magnetic-field configurations that can be stretched without being destroyed 
by the concurrent refinement of the reversal scale and that add up to give rise 
to exponential growth of the magnetic energy. Below we give an analytical 
demonstration of this. It is a (somewhat modified) version of an ingenious ar- 
gument originally proposed by Zeldovich et al. in 1984 |71| . A reader looking 
solely for a broad qualitative picture of isotropic MHD turbulence may skip 
to S IO 

Overcoming diffusion. Let us study magnetic fields with reversals at subviscous 
scales: at these scales, the velocity field is smooth and can, therefore, be expanded 

(f , x) = u' [t, 0) + o'^{t)x'^ + . . . , (37) 

where (T^(i) is the rate-of-strain tensor. The expansion is around some reference 
point X = 0. We can always go to the reference frame that moves with the velocity 

This is not the only existing way of diagnosing the field structure. Ott and cowork- 
ers studied field reversals by measuring magnetic-flux cancellations [JS] . Chertkov 
et al. [7] considered two-point correlation functions of the magnetic field in a 
model of small-scale dynamo and found large-scale correlations along the field 
and short-scale correlations across. 
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at this point, so that u^{t,0) = 0. Let us seek the solution to Eq. ^ with velocity 
H37p as a sum of random plane waves with time-dependent wave vectors: 



B'(t,x) 



d^ko 
(27r)3 



B'{t,ko)e 



ik(t,ko ) 'X 



(38) 



where k(0, ko) = ko, so -B'(0,ko) — i3o(ko) is the Fourier transform of the initial 
field. Since Eq. ([2]) is linear, it is sufficient to ensure that each of the plane waves 
is individually a solution. This leads to two ordinary differential equations for every 
ko: 



dtB' 



dtki 



(39) 



subject to initial conditions -B'(0,ko) = -Bo(ko) and fci(0,ko) = koi. The solution of 
these equations can be written explicitly in terms of the Lagragian transformation 
of variables xo x(t,xo), where 



9ta;'(t,xo) = u{t,x{t,xo)) = (T^(t)a::"'(t, xo), a::'(0,xo) = x},. 



(40) 



Because of the linearity of the velocity field, the strain tensor dx^/dxo^ and its 
inverse dx^/dx'' are functions of time only. At t = 0, they are unit matrices. At 
t > 0, they satisfy 



dt 



dx' 
dx^ 



dxv;' 



Q 

dx^ 



_,dx[_ 
dxi 



We can check by direct substitution that 



dx' 
dxV^ 



Bo"'(ko)exp 



T) I dt'P{t') 



fci(t, ko) 



dxp 
dx^ 



kor 



(41) 



(42) 



These formulae express the evolution of one mode in the integral (|38p . Using the 
fact that det{dxQ/dx') = 1 in an incompressible flow and that Eq. (|42|l therefore 
establishes a one-to-one correspondence k <-> ko, it is easy to prove that the volume- 
integrated magnetic energy is the sum of the energies of individual modes: 



{B')it) 



d^x\B{t,x)f 



d^ko 

(2^ 



B(t,ko)r 



From Eq. (ggj. 



|B(t,ko)!^ = Bo(ko) • M(t) ■ B,*(ko) exp 



-2r, 



f dt'ku 
Jo 



where the matrices M and M ^ have elements defined by 



dx^ dx' 
dxfdx^ 



and M""(f) = 



Oxq dxQ 
dx^ dx^ ' 



(43) 



(44) 



(45) 



respectively. They are the co- and contravariant metric tensors of the inverse La- 
grangian transformation x — > xq. 

Let us consider the simplest possible case of a flow (|37p with constant a = 
diag {Ai, A2, A3}, where Ai > A2 > > A3 and Ai -|- A2 -|- A3 = by incompressibility. 
Then M = diag {e^^i*, e^^^*, e^^^tj ^nd Eq. dH} becomes, in the limit t 00, 
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|B(t,ko)l' ~ |B,i(ko)l"exp 



Ai A2 A3 



(46) 



where we have dropped terms that decay exponentially with time compared to those 
retainedF^ We see that for most ko, the corresponding modes decay superexponen- 
tially fast with time. The domain in the ko space containing modes that are not 
exponentially small at any given time t is given by 



\lt/r^ XiX2t/rj Xi\X3\te^^^yrj 

The volume of this domain at time t is ~ Xl(X2\X:i\y^^{t/r]y^^^e^'^K Within this 
volume, |B(t,ko)|^ ~ |Bo(ko)pe^^i*. Using A3 = -Ai - A2 and Eq. we get 

(B')(t)ocexp[(Ai-A2)t]. (48) 

Let us discuss the physics behind the Zeldovich et al. calculation sketched above. 
When the magnetic field is stretched by the flow, it naturally aligns with the stretch- 
ing Lyapunov direction: B ~ eiBge'^i*. The wave vector k has a tendency to align 
with the compression direction: k ~ e3fco3e''^^'*, which makes most modes decay 
superexponentially. The only ones that survive are those whose ko's were nearly 
perpendicular to 63, with the permitted angular deviation from 90° decaying ex- 
ponentially in time ~ e~'^^'*. Since the magnetic field is solenoidal, Bo -L ko, the 
modes that get stretched the most have Bo || ei and ko || 62 (Fig. [5]). In contrast, 
in 2D, the field aligns with ei and must, therefore, reverse along 62, which is always 
the compression direction (Fig. [5} , so the stretching is always overwhelmed by the 
diffusion and no dynamo is possible (as should be the case according to the rigourous 
early result of Zeldovich [7Up. 

The above construction can be generalised to time-dependent and random ve- 
locity fields. The matrix M is symmetric and can, therefore, be diagonalised by 
an appropriate rotation R of the coordinate system: M = R"^ • L ■ R, where, by 
definition, L = diag |e^i'*\ e^^''\ e^^'*' | . It is possible to prove that, as t ^ 00, 

R(t) {61,62,63} and C,i{t)/2t Xi, where are constant orthogonal unit vec- 
tors, which make up the Lyapunov basis, and Ai are the Lyapunov exponents of 
the flow [21]. The instantaneous values of C,i{t)/2t are called finite-time Lyapunov 
exponents. For a random flow, C,i{t) axe random functions. Eq. (|48|) generalises to 



(B2)(f)cxexp[(Ci-C2)/2], (49) 

where the overline means averaging over the distribution of C,i. The only random 
flow for which this distribution is known is a Gaussian white-in-time velocity first 
considered in the dynamo context in 1967 by Kazantsev [311 The distribution of 
C,i for this flow is Gaussian in the long-time limit and Eq. (|49|l gives (B^) oc g'^'^^-'^^*. 



If A2 — 0, ko2/X2 in Eq. (|46p is replaced with 2^02^- The case A2 < is treated in 
a way similar to that described below and also leads to magnetic-energy growth. 
The difference is that for A2 > 0, the magnetic structures are flux sheets, or 
ribbons, while for A2 < 0, they are flux ropes. 

As the only analytically solvable model of random advection, Kazantsev's model 
has played a crucial role. Developed extensively in 1980s by Zeldovich and cowork- 
ers [72], the model became a tool of choice in the theories of anomalous scaling 
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Fig. 5. Magnetic fields vs. the Lyapunov directions (from [57] )• Zeldovich et al. 
|71) did not give this exact interpretation of their calculation because the folded 
structure of the field was not yet clearly understood at the time. 



where Ai — {(^i)/2t [7j. For Kazantsev's velocity, it is also possible to calculate the 
magnetic-energy spectrum [311 136) , which is the spectrum of the direction reversals. 
It has a peak at the resistive scale and a k^^^^ power law stretching across the sub- 
viscous range, <^ k <^ . This scaling appears to be corroborated by numerical 
simulations [571127] . 

Folded structure revisited. We shall now give a very simple demonstration that 
linear stretching does indeed produce folded fields with straight/curved field lines 
corresponding to larger/smaller field strength. Using Eq. ([2]), we can write evolution 
equations for the field strength B = |B|, the field direction b = B/_B and the 
field-line curvature K = b • Vb. Omitting the resistive terms. 



dB_ 
H 
dh 

m 

dK 

~dt 



= (bb:Vu)B, (50) 

= b- (Vu) ■ (I -bb), (51) 

= K ■ (Vu) ■ (I - bb) - 2(bb : Vu)K - [b • (Vu) • K] b 

+bb : (VVu) ■ (I -bb). (52) 



For simplicity, we again use the velocity field (|37|l with constant o — diag {Ai, A2, A3}. 
Then the stable fixed point of Eq. (|51|l in the comoving frame is b = ei (magnetic 
field aligns with the principal stretching direction), whence B oc e^^*. Since K-b — 0, 
we set Ki = 0. Denoting S — bb : Wu, we can now write Eq. (|52|) as 



dK2 

dt 



dt 



(3Ai+A2)A'3 + X'3 



(53) 



and intermittency that flourished in 1990s [16] (in this context, it has been associ- 
ated with the name of Kraichnan who, independently from Kazantsev, proposed 
to use it for the passive scalar problem [35] )• It remains useful to this day as old 
theories are reevaluated and new questions demand analytical answers |55j . 
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G. K. Batdielor (1920-2000) A. Schliiter L. Biermann (1907-1986) 

Fig. 6. Photo of A. Schliiter courtesy of MPI fiir Plasmaphysik; photo of L. Bier- 
mann courtesy of Max-Planck-Gesellschaft/AIP Emilio Segre Visual Archives. 

Both components of K decay exponentially^ until they are comparable to the 
inverse scale of the velocity field (i.e., the terms containing Wu become important). 
The stationary solution is K2 ~ X'2/(2Ai — A2), — X'3/(3Ai -|- A2). 

If the field is to reverse direction, it must turn somewhere (see Fig. [S}. At such 
a turning point, the field must be perpendicular to the stretching direction. Setting 
61 = 0, we find two fixed points of Eq. (|51|l under this condition: b = 62 and b = 63. 
Only the former is stable, so the field at the turning point will tend to align with the 
"nuU" direction. Thus, stretching favours configurations with field reversals along 
the "null" direction, which are also those that survive diffusion (see Fig. [5]). From 
Eq. (|52p we find that at the turning point, K2 = 0, K3 = X'3/(Ai + 3A2), while Ki 
grows at the rate Ai — 2A2 (assumed positive). This growth continues until limited 
by diffusion at AT ~ l/^>7- The strength of the field in this curved region is oc e^^*, 
so the fields are weaker than in the straight segments, where B oc e"^^'. 

When the problem is solved for the Kazantsev velocity, the above solution gener- 
alises to a field of random curvatures anticorrelated with the magnetic-field strength 
and with a stationary PDF of K that has a peak at K ^ flow scale"^ and a power 
tail ~ K~^^^^ describing the distribution of curvatures at the turning points [55j . 
Numerical simulations support these results [57) . 



3.2 Saturation of the dynamo 

The small-scale dynamo gave us exponentially growing magnetic fields with 
energy concentrated at small (resistive) scales. How is the growth of magnetic 
energy saturated and what is the final state? Will magnetic energy stay at 
small scales or will it proceed to scale-by-scale equipartition via some form 
of inverse cascade? This basic dichotomy dates back to the 1950 papers by 
Batchelor [2] and Schliiter & Biermann [6D] . Batchclor thought that magnetic 
field was basically analogous to the vorticity field a; = V x u (which satisfies 
the same equation ^ except for the difference between 1] and v) and would, 

K3 decays faster than K2. If the velocity is exactly linear {S = 0), K aligns 
with 62 and decreases indefinitely, while the combination 1/(2^-^2/^1) g^^yg 
constant. This rhymes with the result that can be proven for a linear Kazantsev 
velocity: at zero 77, B oc e^i/^, while BK^^^ oc .^j^^ X2 = {(,2) /2t = 0. 
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therefore, saturate at a low energy. (B^) ^ i?e^^/^(u^), with a spectrum 
peaked at the viscous scale. Schliiter and Biermann disagreed and argued that 
the saturated state would be a scale-by-scale balance between the Lorentz and 
inertial forces, with turbulent motions at each scale giving rise to magnetic 
fluctuations of matching energy at the same scale. Schliiter and Biermann's 
argument (and, implicitly, also Batchelor's) was based on the assumption of 
locality of interaction (in scale space) between the magnetic and velocity 
fields: locality both of the dynamo action and of the back reaction. As we 
saw in the previous section, this assumption is certainly incorrect for the 
dynamo: a linear velocity field, i.e., a velocity field of a formally infinitely large 
(in practice, viscous) scale can produce magnetic fields with reversals at the 
smallest scale allowed by diffusion. A key implication of the folded structure 
of these fields concerns the Lorentz force, the essential part of which, in the 
case of incompressible flow, is the curvature force B^K.. Since it is a quantity 
that depends only on the parallel gradient of the magnetic field and does not 
know about direction reversals, it will possess a degree of velocity-scale spatial 
coherence necessary to oppose stretching. Thus, a field that is formally at the 
resistive scale will exert a back reaction at the scale of the velocity field. In 
other words, interactions are nonlocal: a random fiow at a given scale 
having amplified the magnetic fields at the resistive scale Ijj ^ I, will see 
these magnetic fields back react at the scale I. Given the nonlocality of back 
reaction, we can update Batchelor's and Schliiter & Biermann's scenarios for 
saturation in the following way [56l [57] . 

The magnetic energy is amplified by the viscous-scale motions until the 
field is strong enough to resist stretching, i.e., until B- VB ~ u- Vu ^ Suf /l,^. 
Since B • VB - B'^K - B'^ /ly (folded field), this happens when 

{B^) ^ 5ul^ Re-^'^iu"). (54) 

Let us suppose that the viscous motions arc suppressed by the back reaction, 
at least in their ability to amplify the field. Then the motions at larger scales in 
the inertial range come into play: while their rates of strain and, therefore, the 
associated stretching rates are smaller than that of the viscous-scale motions, 
they are more energetic [see Eq. Q], so the magnetic field is too weak to resist 
being stretched by them. As the field continues to grow, it will suppress the 
motions at ever larger scales. If we define a stretching scale ls{t) as the scale 
of the motions whose energy is 5ui^ ~ (i?^)(t), we can estimate 

|(B^)^^(i?2)^^~e = const {B^m^et. (55) 

Thus, exponential growtli gives way to secular growth of the magnetic en- 
ergy. This is accompanied by elongation of the folds (their length is always of 
the order of the stretching scale, Z|| ~ Is), while the resistive (reversal) scale 
increases because the stretching rate goes down: 

Z|| (t) ^ Ut) ^ Sul/e ^ Vii3/2^ _ [rj/{5uijk)f^ ^ ^t. (56) 
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This secular stage can continue until the entire incrtial range is suppressed, 
Is ~ i, at which point saturation must occur. This happens after t 
(T^l'^lp-I'^ ^ L/6ul. Using Eqs. we have, in saturation, 

{B^)^{u^), k^L, l^^{^l[&u^lL)\^l^ ^R-^I-^L. (57) 

Comparing Eqs. ((57)) and ((36)l . we see that the resistive scale has increased only 
by a factor of i?e^^'* over its value in the weak-field growth stage. Note that this 

imposes a very stringent requirement on any numerical experiment striving 

^ 1 
to distinguish between the viscous and resistive scales: Pr^a ^ Re. ' ^1. 

If, as in the above scenario, the magnetic field retains its folded struc- 
ture in saturation, with direction reversals at the resistive scale, this explains 
qualitatively why the numerical simulations of the developed isotropic MHD 
turbulence with Pr^ > 1 |57] show the magnetic-energy pile-up at the small 
scales. What then is the saturated state of the turbulent velocity field? We as- 
sumed above that the incrtial-range motions were "suppressed" — this applied 
to their ability to amplify magnetic field, but needed not imply a complete 
evacuation of the inertial range. Indeed, simulations at modest Pr™ show a 
powerlike velocity spectrum [S7l [57] . The most obvious class of motions that 
can populate the inertial range without affecting the magnetic-field strength 
are a type of Alfven waves that propagate not along a mean (or large-scale) 
magnetic field but along the folded structure (Fig. [7]). Mathematically, the 
dispersion relation for such waves is derived via a linear theory carried out 
for the inertial-range perturbations (L^^ <C fc <C 1^^) of the tensor BiBj (cf. 
[55]). The unperturbed state is the average of this tensor over the subviscous 
scales: {BiBj) = bibj{B^), where (B^) is the total magnetic energy and the 
tensor bibj only varies at the outer scale L [Eq. (|57|) ]. The resulting disper- 
sion relation is w = ±|k- b|(i3^)^/^ [SSj. The presence of these waves will not 
change the rcsistive-scale-dominatcd nature of the magnetic-energy spectrum, 
but should be manifest in the kinetic-energy spectrum. There is, at present, 
no theory of a cascade of such waves, although a line of argument similar to §[2] 
might work, since it docs not depend on the field having a specific direction. 
A numerical detection of these waves is also a challenge for the future. 

What we have proposed above can be thought of as a modernised ver- 
sion of the Schliiter & Biermann scenario, retaining the intermediate secular- 
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growth stage and saturation with {B'^) ^ {u'^), but not scalc-by-scale equipar- 
tition. However, an alternative possibihty, which is in a similar relationship 
to Batchclor's scenario, can also be envisioned. In Eqs. ([55]) - ([551) . the scale 
at which diffusion cuts off the small-scale magnetic fluctuations was assumed 
to be determined by the stretching rate Sui^/lg. However, since the nonlinear 
suppression of the viscous-scale eddies only needs to eliminate motions with 
bb : Vu 7^ [Eq. ([SO]) ]. two-dimensional "interchange" motions (velocity gra- 
dients _L b) are, in principle, allowed to survive at the viscous scale. These 
could "two-dimensionally" mix the direction-reversing magnetic fields at the 
rate Sui^Jlu — much faster tlian tlie unsuppressed larger-scale stretching can 
amplify the field, — with tlie consequence that the resistive scale is pinned 
at the value given by Eq. ([55)1 and the field cannot grow above the Batch- 
elor limit (|54[) . The mixing efficiency of the suppressed motions is the key 
to choosing between the two saturation scenarios. Numerical simulations [57] 
corroborate the existence of an intermediate stage of slower-than-exponential 
growth accompanied by fold elongation and a modest increase of the resistive 
scale [Eq. This tips the scales in favour of the first scenario, but, in view 

of limited resolutions, we hesitate to declare the matter definitively resolved. 

3.3 Turbulence and magnetic fields in galaxy clusters 

The intracluster medium (hereafter, ICM) is a hot (T ~ 10^ K) diffuse 
{n ^ 10^^. ..10^'^ cm^"^) fully ionised plasma, which accounts for most of the 
luminous matter in the Universe (note that it is not entirely dissimilar from 
the ionised phases of the interstellar medium: e.g., the so-called hot ISM). 
It is a natural astrophysical environment to which the large-Prm isotropic 
regime of MHD turbulence appears to be applicable: indeed, Eq. ([55)) gives 
Prm ~ 1029. 

The ICM is believed to be in a state of turbulence driven by a variety 
of mechanisms: merger events, galactic and subcluster wakes, active galac- 
tic nuclei. One expects the outer scale L ^ kpc and the velocity 
dispersions 5ul ^ 10^... 10"^ km/s (a fraction of the sound speed). Indirect 
observational evidence supporting the possibility of a turbulent ICM with 
roughly these parameters already exists (an apparently powcrlike spectrum 
of pressure fluctuations found in the Coma cluster [61] , broadened abundance 
profiles in Perseus believed to be caused by turbulent diffusion [53]), and di- 
rect detection may be achieved in the near future [29] . However, there is as 
yet no consensus on whether turbulence, at least in the usual hydrodynamic 
sense, is a generic feature of clusters [15]. The main difficulty is the very large 
values of the ICM viscosity obtained via the standard estimate v ~ i'th,i'^mfp, 
where Vth,i ~ 10'^ km/s is the ion thermal speed and A^fp ^ 1...10 kpc is the 
ion mean free path. This gives Re ~ 10^ if not less, whicli makes the exis- 
tence of a well-developed inertial range problematic. Postponing the problem 
of viscosity until §|4l we observe that the small-scale dynamo does not, in 
fact, require a turbulent velocity field in the sense of a broad inertial range: 
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in the weak-field regime discussed in §13111 the dynamo was controUed by the 
smooth single-scale random flow associated with the viscous-scale motions; 
in saturation, we argued in § 13.21 that the main effect was the direct nonlocal 
interaction between the outer-scale (random) motions and the magnetic flcld. 
Given the available menu of large-scale stirring mechanisms in clusters, it is 
likely that, whatever the value of Re, the velocity field is randomF^ 

The cluster turbulence is certainly magnetic. The presence of magnetic 
fields was first demonstrated for the Coma cluster, for which Willson detected 
in 1970 a diffuse synchrotron radio emission [55] and Kim et al. in 1990 were 
able to estimate directly the magnetic-field strength and scale using the Fara- 
day rotation measure (RM) data [32] ■ Such observations of magnetic fields in 
clusters have now become a vibrant area of astronomy (reviewed most recently 
in [25]), usually reporting a field B ~ 1...10 at scales ~ 1...10 kpc [TUjllj 
All of this field is small-scale fluctuations: no appreciable mean component 
has been detected. The field is dynamically significant: the magnetic energy 
is less but not much less than the kinetic energy of the turbulent motions. 

Do clusters fit the theoretical expectations reviewed above? The magnetic- 
field scale seen in clusters is usually 10 to 100 times smaller than the expected 
outer scale of turbulent motions and, indeed, is also smaller than the viscous 
scale based on Re ^ 10^. However, it is certainly far above the resistive scale, 
which turns out be ~ 10'^. ..10^ km! Faced with these numbers, we must 
suspend the discussion of cluster physics and finally take account of the fact 
that astrophysical bodies are made of plasma, not of an MHD fiuid. 

4 Enter plasma physics 
4.1 Braginskii viscosity 

In all of the above, we have used the MHD equations H]) and ([2]) to develop 
turbulence theories supposed to be relevant for astrophysical plasmas. His- 
torically, such has been the approach followed in most of the astrophysical 
literature. The philosophy underpinning this approach is again that of univer- 
sality: the "microphysics" at and below the dissipation scale are not expected 
to matter for the fiuid-like dynamics at larger scales. However, in considering 
the MHD turbulence with large Pr^, we saw that dissipation scales, deter- 
mined by the values of the viscosity v and magnetic diffusivity 77, played a 
very prominent role: the growth of the small-scale magnetic fields was con- 
trolled by the turbulent rate of strain at the viscous scale and resulted in the 
magnetic energy piling up, in the form of direction-reversing folded fields, at 

Numerical simulations of the large- Pr^ regime at currently accessible resolutions 
also rely on a random forcing to produce "turbulence" with Re ~ 1...10'^ [57ll27j . 
^® Because of the availability of the RM maps from extended radio sources in clus- 
ters, it is possible to go beyond field-strength and scale estimates and construct 
magnetic-energy spectra with spatial resolution of ~ 0.1 kpc |68) . 
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the resistive scale — both in the growth and saturation stages of the dynamo. 
It is then natural to revisit the question of whether the Laplacian diffusion 
terms in Eqs. HI)-© arc a good description of the dissipation in astrophysical 
plasmas. 

The answer to this question is, of course, that they are not. A necessary 
assumption in the derivation of these terms is that the ion cyclotron frequency 
f2i = eB/rriiC exceeds the ion-ion collision frequency vu or, cquivalcntly, the 
ion gyroradius pi ~ vth.i/^^i exceeds the mean free path Amfp = vth,i/i^ii- This 
is patently not the case in many astrophysical plasmas: for example, in galaxy 
clusters, Amfp ~ 1...10 kpc, while pi ~ 10^ km. In such a weakly collisional 
magnetised plasma, the momentum equation ([T|) assumes the following form, 
valid at spatial scales ^ pi and at time scales ^ ^i^^, 

^ = -^{p^ + ^)+^- [bb(Pi - P\\ + B')] + f , (58) 

where p± and p\\ are plasma pressures perpendicular and parallel to the local 
direction of the magnetic field, respectively, and we have used B • VB = 
V • (hhB^). The evolution of the magnetic field is controlled by the electrons 
— the field remains frozen into the flow and we may use Eq. ^ with 77 = 0. 

If we arc interested in subsonic motions, V(px + B"^ /2) in Eq. ([55)1 can be 
found from the incompressibility condition V • u = and the only quantity 
still to be determined is p± —p\\ - The proper way to compute it is by a rather 
lengthy kinetic calculation due to Braginskii [5], which cannot be repeated 
here. The result of this calculation can, however, be obtained in the following 
heuristic way [55] . 

The fundamental property of charged particles moving in a magnetic field 
is the conservation of the first adiabatic invariant p ~ miV^/2B^]\ When 
Amfp S> Pi, this conservation is only weakly broken by collisions. As long as p 
is conserved, any change in B must be accompanied by a proportional change 
in p±. Thus, the emergence of the pressure anisotropy is a natural consequence 
of the changes in the magnetic-field strength and vice versa: indeed, summing 
up the first adiabatic invariants of all particles, we get p±/B = const. Then 

1 dpj_ 1 dB P^-P\\ 

i^ii (59) 



p± dt B dt pj 

where the second term on the right-hand sight represents the collisional relax- 
ation of the pressure anisotropy p±^ — p\\ at the rate ^ Wth.i/ Amfp 111 Using 
Eq. ([50)) for B and balancing the terms in the rhs of Eq. ([50]) , we get 



It may be helpful to the reader to think of this property as the conservation of 
the angular momentum of a gyrating particle: miV±pi oc miv\/B = 2p. 
This is only valid if the characteristic parallel scales of all fields are larger than 
Amfp. In the colhsionless regime, fcyAmfp ^ 1, we may assume that the pressure 
anisotropy is relaxed in the time particles streaming along the field cover the 
distance k7^: this entails replacing va in Eq. (|59p by fcyUth.i. 
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P± - P\\ = J^ll = '^llbb : Vu, (60) 

where i^y p/i^u ~ Vth,i^niip is the "parallel viscosity." This equation turns 
out to be exact [5] up to numerical prefactors in the definition of i/y . 
The energy conservation law based on Eqs. ([55)1 and ^ is 

where Ohmic diffusion has been omitted. Thus, the Braginskii viscosity only 
dissipates such velocity gradients that change the strength of the magnetic 
field. The motions that do not affect B are allowed to exist in the subviscous 
scale range. In the weak-field regime, these motions take the form of plasma 
instabilities. When the magnetic field is strong, a cascade of shcar-Alfvcn 
waves can be set up below the viscous scale. Let us elaborate. 

4.2 Plasma instabilities 

The simplest way to see that the pressure anisotropy in Eq. ((58|) leads to 
instabilities is as follows [58]. Imagine a "fluid" solution with u. p±, B 
changing on viscous time and spatial scales, t ^ | Vu|^^ ~ li^/Sui^ and I ^ l^. 
Would such a solution be stable with respect to fast {lu ^ | Vu|^-^) small-scale 
(fc ;» perturbations? Linearising Eq. ([55]) and denoting perturbations by 
S, we get 

- iuSu = -ik {Spi_ + B5B) + {p±_ - p|| + B^) (5K 

+ ih fc|| [(5p_L - 5p\\ - (p± - P|| - B^) 5B/B] , (62) 

where the perturbation of the field curvature is (5K = fc^(5u^/iw [see Eq. ([5^ ]. 
We see that regardless of the origin of the pressure anisotropy, the shear- 
Alfven-polarised perturbations {5\i oc k x b) have the dispersion relation 

w = ±fc|| (pi +5^)1/2^ (63) 

When p|| — > B^, w is purely imaginary and we have what is known as 
the firehose instability (HH [SI [^ [17] . The growth rate of the instability is 
(X k\\ , which means that the fastest-growing perturbations will be at scales far 
below the viscous scale or, indeed, the mean free path. Therefore, adopting 
the Braginskii viscosity [Eq. ([60]) ] exposes a fundamental problem with the 
use of the MHD approximation for fully ionised plasmas: the equations are 
ill posed wherever p\\ — p± > B^. To take into account the instability and its 
impact on the large-scale dynamics, the fluid equations must be abandoned 
and a kinetic description adopted. A linear kinetic calculation shows that the 
instability growth rate peaks at k^^pi ^ 1, so the fluctuations grow fastest at 
the ion gyroscale. While the flrehose instability occurs in regions where the 
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velocity field leads to a decrease in the magnetic- field strength [Eq. (pD)) ]. a 
kinetic calculation of the pressure perturbations in Eq. (|62p shows that another 
instability, called the mirror mode |51| , is triggered wherever the field increases 
(p_L > Its growth rate is also oc fcy and peaks at the ion gyroscale. 

In weakly collisional astrophysical plasmas such as the ICM, the random 
motions produced by the large-scale stirring will stretch and fold magnetic 
fields, giving rise to regions both of increasing and decreasing field strength 
(^ I3.ip . The instabilities should, therefore, be present in weak-field regions 
where — p±\ > and, since their growth rates are much larger than 
the fluid rates of strain, their growth and saturation should have a profound 
effect on the structure of the turbulence. A quantitative theory of what ex- 
actly happens is not as yet available, but one might plausibly expect that the 
fluctuations excited by the instabilities will lead to some effective renormali- 
sation of both the viscosity and the magnetic diffusivity. A successful theory 
of turbulence in clusters requires a quantitative calculation of this effective 
transport. In particular, this should resolve the uncertainties around the ICM 
viscosity and produce a prediction of the magnetic-fleld scale to be compared 
with the observed values reviewed in S I3.3I 

In the solar wind, the plasma is magnetised {pi ^ 10^ km), while collisions 
are virtually absent: the mean free path exceeds the distance from the Sun 
(10^ km). Ion pressure (temperature) anisotropics with respect to the field 
direction were directly measured in 1970s [T7l[4T]. As was first suggested by 
Parker [51] , firehose and mirror instabilities (as well as several others) should 
play a major, although not entirely understood role [19] . A vast geophysical 
literature now exists on this subject, which cannot be reviewed here. 

4.3 Kinetic turbulence 

The instabilities are quenched when the magnetic field is sufficiently strong: 
overwhelms in the second term on the right-hand side of Eq. (j58p . 

If we use the collisional estimate (|60|) . this happens when 

B'^ ^ v\\5uiJK Re~^''^5ul. (64) 

The firehose- unstable perturbations become Alfvcn waves in this limit. In 
the strong-field regime {6B <C i?o), the appropriate mathematical description 
of the weakly collisional turbulence of Alfven waves is the low-frequency limit 
of the plasma kinetic theory called the gyrokinetics [TSl [53] It is obtained 
under an ordering scheme that stipulates 

k^p^^l, uj/f2i ^ k\i/kj_ ^ Su/vA ^ SB/Bq <t: 1. (65) 

^® While originally developed and widely used for fusion plasmas, this "kinetic-fiuid" 
description has only recently started to be applied to astrophysical problems such 
as the relative heating of ions and electrons by Alfvenic turbulence in advection- 
dominated accretion flows [52j . 
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The second relation in Eq. coincides with the GS critical-balance conjec- 
ture (|14p if the latter is treated as an ordering assumption. The gyrokinetics 
can be cast as a systematic expansion of the full kinetic description of the 
plasma in the small parameter e ~ fc||/fc± — the direct generalisation of the 
similar expansion of MHD equations given at the end of ^ \2A[ It turns out 
that the decoupling of the Alfven-wave cascade that we demonstrated there is 
also a property of the gyrokinetics and that this cascade is correctly described 
by the RMHD equations ([22|) and ([23|l all the way down to the ion gyroscale, 
k±pi ~ 1 [S5]. Broad fluctuation spectra observed in the solar wind [H] and 
in the ISM [1] are likely to be manifestations of just such a cascade. The slow 
waves and the entropy mode are passively mixed by the Alfven-wave cascade, 
but Eqs. ((29|) - ([3T|) have to be replaced by a kinetic equation. 

When magnetic fields are not stronger than the turbulent motions — as is 
the case for clusters, where the magnetic energy is, in fact, quite close to the 
threshold (|64p — the situation is more complicated and much more obscure 
because small-scale dynamo (ii l3.1|) . back reaction (§[3?2]), plasma instabilities 
(in weak-field regions such as, for example, the bending regions of the folded 
fields, i) l3.ip . and Alfven waves (possibly of the kind discussed in ^ 13. 2p all 
enter into the mix and remain to be sorted out. 

5 Conclusion 

We conclude here in the hope that we have provided the reader with a fair 
overview of the state of affairs to which the MHD turbulence theory has ar- 
rived after its first fifty years. Perhaps, despite much insight gained along the 
way, not very far. It is clear that a simple extension of Kolmogorov's theory 
has so far proven unattainable. Two of the key assumptions of that theory 
— isotropy and locality of interactions — are manifestly incorrect for 
MHD. Indeed, even the applicability of the fundamental principle of small- 
scale universality is suspect. Although a fair amount is known about Alfvenic 
turbulence, progress in answering many important astrophysical questions (see 
the Introduction) has been elusive because there is little knowledge of the gen- 
eral spectral and structural properties of the fully developed turbulence in an 
MHD fluid and, more generally, in magnetised weakly coUisional or collision- 
less plasmas. Thus, while many unanswered questions demand further effort 
on MHD turbulence, there is also an imperative, mandated by astrophysical 
applications, to go beyond the fluid description. 

Acknowledgments. We would like to thank Russell Kulsrud who is origi- 
nally responsible for our interest in these matters. We are grateful to N. A. 
Lipunova, M. Bernard, R. A. Sunyaev, R. Buchstab and H. Lindsay for help- 
ing us obtain the photos used in this chapter. Our work was supported by 
the UKAFF Fellowship, PPARC Advanced Fellowship, King's College, Cam- 
bridge (A.A.S.) and in part by the US DOE Center for Multiscale Plasma 
Dynamics (S.C.C.). 



Turbulence and Magnetic Fields in Astrophysical Plasmas 



27 



References 

1. Armstrong JW, Rickett BJ, Spangler SR (1995) Electron density power spec- 
trum in the local interstellar medium. Astrophys J 443:209-221 

2. Batchelor GK (1950) On the spontaneous magnetic field in a conducting liquid 
in turbulent motion. Proc Roy Soc London A 201:405-416 

3. Belcher JW, Davis L (1971) Large-amplitude Alfven waves in interplanetary 
medium, 2. J Geophys Res 76:3534-3563 

4. Biskamp D, Miiller W-C (2000) Scaling properties of three-dimensional 
isotropic magnetohydrodynamic turbulence. Phys Plasmas 7:4889-4900 

5. Braginskii SI (1965) Transport processes in a plasma. Rev Plasma Phys 1:205- 
310 

6. Chandrasekhar S, Kaufman AN, Watson KM (1958) The stability of the pinch. 
Proc R Soc London A 245:435-455 

7. Chertkov M, Falkovich G, Kolokolov I, Vergassola M (1999) Small-scale turbu- 
lent dynamo. Phys Rev Lett 83:4065-4068 

8. Childress S, Gilbert AD (1995) Stretch, Twist, Fold: The Fast Dynamo. 
Springer 

9. Cho J, Lazarian A, Vishniac ET (2002) Simulations of magnetohydrodynamic 
turbulence in a strongly magnetized medium. Astrophys J 564:291-301 

10. Clarke TE, Kronberg PP, Bohringer H (2001) A new radio-X-ray probe of 
galaxy cluster magnetic fields. Astrophys J 547:L111-L114 

11. Coleman PJ (1968) Turbulence, viscosity, and dissipation in the solar-wind 
plasma. Astrophys J 153:371-388 

12. Dennis TJ, Chandran BDG (2005) Turbulent heating of galaxy-cluster plasmas. 
Astrophys J 622:205-216 

13. Dobrowolny M, Mangeney A, Veltri P (1980) Fully developed anisotropic hy- 
dromagnetic turbulence in the interplanetary space. Phys Rev Lett 45:144-147 

14. Elsasser WM (1950) The hydromagnetic equations. Phys Rev 79:183 

15. Fabian AC, Sanders JS, Crawford CS, Conselice CJ, Gallagher JS, Wyse RFC 
(2003) The relationship between optical Ha filaments and the X-ray emission 
in the core of the Perseus cluster. Mon Not R Astr Soc 344:L48-L52 

16. Falkovich G, Gaw§dzki K, Vergassola M (2001) Particles and fields in fluid 
turbulence. Rev Mod Phys 73:913-975 

17. Feldman WC, Asbridge JR, Bame SJ, Montgomery MD (1974) Interpenetrating 
solar wind streams. Rev Geophys Space Phys 12:715-723 

18. Frieman EA, Chen L (1982) Nonlinear gyrokinetic equations for low-frequency 
electromagnetic waves in general plasma equilibria. Phys Fluids 25:502-508 

19. Gary SP (1993) Theory of space plasma microinstabilities. Cambridge: Cam- 
bridge University Press 

20. Galtier S, Nazarenko SV, Newell AC, Pouquet A (2000) A weak turbulence 
theory for incompressible magnetohydrodynamics. J Plasma Phys 63:447-488 

21. Goldhirsch I, Sulem P-L, Orszag SA (1987) Stability and Lyapunov stability of 
dynamical systems: a differential approach and a numerical method. Physica D 
27:311-337 

22. Goldreich P, Sridhar S (1995) Toward a theory of interstellar turbulence. II. 
Strong Alfvenic turbulence. Astrophys J 438:763-775 

23. Goldreich P, Sridhar S (1997) Magnetohydrodynamic turbulence revisited. As- 
trophys J 485:680-688 



28 Alexander A Schekochihin and Steven C Cowley 

24. Goldstein ML, Roberts DA (1999) Magnetohydrodynamic turbulence in the 
solar wind. Phys Plasmas 6:4154-4160 

25. Govoni F, Feretti L (2004) Magnetic fields in clusters of galaxies. Int J Mod 
Phys D 13:1549-1594 

26. Gruzinov AV, Diamond PH (1996) Nonlinear mean field electrodynamics of 
turbulent dynamos. Phys Plasmas 3:1853-1857 

27. Haugen NEL, Brandenburg A, Dobler W (2004) Simulations of nonhelical hy- 
dromagnetic turbulence. Phys Rev E 70:016308 

28. Higdon JC (1984) Density fluctuations in the interstellar medium: evidence for 
anisotropic magnetogasdynamic turbulence. I. Model and astrophysical sites. 
Astrophys J 285:109-123 

29. Inogamov NA, Sunyaev RA (2003) Turbulence in clusters of galaxies and X- 
ray-line profiles. Astron Lett 29:791-824 

30. Iroshnikov RS (1963) Turbulence of a conducting fluid in a strong magnetic 
fleld. Astron Zh 40:742-750 [English translation: Sov Astron 7:566-571 (1964) 
Note that Iroshnikov's initials are given incorrectly as PS in the translation] 

31. Kazantsev AP (1967) Enhancement of a magnetic field by a conducting fluid. 
Zh Eksp Teor Fiz 53:1806-1813 [English translation: Sov Phys JETP 26:1031- 
1034 (1968)] 

32. Kim K-T, Kronberg PP, Dewdney PE, Landecker TL (1990) The halo and 
magnetic field of the Coma cluster of galaxies. Astrophys J 355:29-37 

33. Kolmogorov AN (1941) The local structure of turbulence in incompressible 
viscous fluid at very large Reynolds numbers. Dokl Akad Nauk SSSR 30:299- 
303 [English translation: Proc Roy Soc London A 434:9-13 (1991)] 

34. Kraichnan RH (1965) Inertial-range spectrum of hydromagnetic turbulence. 
Phys Fluids 8:1385-1387 

35. Kraichnan RH (1968) Small-scale structure of a scalar fleld convected by tur- 
bulence. Phys Fluids 11:945-953 

36. Kulsrud RM, Anderson SW (1992) The spectrum of random magnetic fields 
in the mean fleld dynamo theory of the galactic magnetic fleld. Astrophys J 
396:606-630 

37. Landau LD, Lifshitz EM (1987) Fluid Mechanics. Butterworth-Heinemann 

38. Lithwick Y, Goldreich P (2001) Compressible magnetohydrodynamic turbu- 
lence in interstellar plasmas. Astrophys J 562:279-296 

39. Lithwick Y, Goldreich P (2003) Imbalanced weak magnetohydrodynamic tur- 
bulence. Astrophys J 582:1220-1240 

40. Maron J, Goldreich P (2001) Simulations of incompressible magnetohydrody- 
namic turbulence. Astrophys J 554:1175-1196 

41. Marsch E, Ao X-Z, Tu C-Y (2004) On the temperature anisotropy of the core 
part of the proton velocity distribution function in the solar wind. J Geophys 
Res 109:A04102 

42. Matthaeus WH, Goldstein ML (1982) Measurement of the rugged invariants of 
magnetohydrodynamic turbulence in the solar wind. J Geophys Res 87:6011- 
1028 

43. Matthaeus WH, Goldstein ML, Roberts DA (1990) Evidence for the presence 
of quasi-two-dimensional nearly incompressible fluctuations in the solar wind. 
J Geophys Res 95:20673-20683 

44. Meneguzzi M, Frisch U, Pouquet A (1981) Helical and nonhelical turbulent 
dynamos. Phys Rev Lett 47:1060-1064 



Turbulence and Magnetic Fields in Astrophysical Plasmas 



29 



45. Montgomery D, Turner L (1981) Anisotropic magnetohydrodynamic turbulence 
in a strong external magnetic field. Phys Fluids 24:825-831 

46. Miiller W-C, Biskamp D, Grappin R (2003) Statistical anisotropy of magnto- 
hydrodynamic turbulence. Phys Rev E 67:066302 

47. Ng CS, Bhattacharjee A (1996) Interaction of shear- Alfven wave packets: im- 
plication for weak magnetohydrodynamic turbulence in astrophysical plasmas. 
Astrophys J 465:845-854 

48. Obukhov AM (1949) Structure of a temperature field in a turbulent fiow. Izv 
Akad Nauk SSSR Ser Geogr Geofiz 13:58-69 

49. Ott E (1998) Chaotic fiows and kinematic magnetic dynamos: A tutorial review. 
Phys Plasmas 5:1636-1646 

50. Parker EN (1958) Dynamics of the interplanetary gas and magnetic fields. As- 
trophys J 128:664-676 

51. Parker EN (1958) Dynamical instability in an anisotropic ionized gas of low 
density. Phys Rev 109:1874-1876 

52. Quataert E, Gruzinov A (1999) Turbulence and particle heating in advection- 
dominated accretion flows. Astrophys J 520:248-255 

53. Rebusco P, Churazov E, Bohringer H, Forman W (2005) Impact of stochas- 
tic gas motions on galaxy cluster abundance profiles. Mon Not R Astr Soc 
359:1041-1048 

54. Rosenbluth MN (1956) Stability of the pinch. Los Alamos Scientific Laboratory 
Report LA-2030 

55. Schekochihin A, Cowley S, Maron J, Malyshkin L (2002) Structure of small- 
scale magnetic fields in the kinematic dynamo theory. Phys Rev E 65:016305 

56. Schekochihin AA, Cowley CS, Hammett GW, Maron JL, McWilliams JC (2002) 
A model of nonlinear evolution and saturation of the turbulent MHD dynamo. 
New J Phys 4:84 

57. Schekochihin AA, Cowley CS, Taylor SF, Maron JL, McWilhams JC (2004) 
Simulations of the small-scale turbulent dynamo. Astrophys J 612:276-307 

58. Schekochihin AA, Cowley CS, Kulsrud RM, Hammett GW, Sharma P (2005) 
Plasma instabilities and magnetic-field growth in clusters of galaxies. Astrophys 
J 629:139-142 

59. Schekochihin AA, Cowley CS, Dorland WD, Hammett GW, Howes GG, 
Quataert E, Tatsuno T (2007) Kinetic and fiuid turbulent cascades in magne- 
tized weakly coUisional astrophysical plasmas. Astrophys J (submitted; E-print 
'arXiv: 0704. 0044) 

60. Schliiter A, Biermann L (1950) Interstellare Magnetfelder. Z Naturforsch 
5a:237-251 

61. Schuecker P, Finoguenov A, Miniati F, Bohringer H, Briel UG (2004) Probing 
turbulence in the Coma galaxy cluster. Astron Astrophys 426:387-397 

62. Shakura NI, Sunyaev RA (1973) Black holes in binary systems. Observational 
appearance. Astron Astrophys 24:337-355 

63. Shebalin JV, Matthaeus WH, Montgomery D (1983) Anisotropy in MHD tur- 
bulence due to a mean magnetic field. J Plasma Phys 29:525-547 

64. Spitzer L (1962) Physics of fully ionized gases. New York: Wiley 

65. Sridhar S, Goldreich P (1994) Toward a theory of interstellar turbulence. I. 
Weak Alfvenic turbulence. Astrophys J 432:612-621 

66. Strauss HR (1976) Nonlinear, three-dimensional magnetohydrodynamics of 
noncircular tokamaks. Phys Fluids 19:134-140 



30 Alexander A Schekochihin and Steven C Cowley 

67. Vedenov AA, Sagdeev RZ (1958) Some properties of the plasma with 
anisotropic distribution of the velocities of ions in the magnetic field. Sov Phys 
Dokl 3:278 

68. Vogt C, Enfilin TA (2005) A Bayesian view on Fraday rotation maps — Seeing 
the magnetic power spectra in galaxy clusters. Astron Astrophys 434:67-76 

69. Willson MAG (1970) Radio observations of the cluster of galaxies in Coma 
Berenices — The 5C4 survey. Mon Not R Astr Soc 151:1-44 

70. Zeldovich YaB (1956) The magnetic field in the two-dimensional motion of a 
conducting turbulent liquid. Zh Exp Teor Fiz 31:154-155 [English translation: 
Sov Phys JETP 4:460-462 (1957)] 

71. Zeldovich YaB, Ruzmaikin AA, Molchanov SA, Sokoloff DD (1984) Kinematic 
dynamo problem in a linear velocity field. J Fluid Mech 144:1-11 

72. Zeldovich YaB, Ruzmaikin AA, Sokoloff DD (1990) The Almighty Chance. 
Singapore: World Scientific 



